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Abstract: We construct a family of non-Gaussian martingales the marginals 
of which are all Gaussian. We give the predictable quadratic variation of 
Q^ these processes and show they do not have continuous paths. These pro- 
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1. Introduction 



This paper is concerned with constructing an alternative process to the Brown- 
ian motion; a martingale with Gaussian marginals, yet not a Gaussian process. 
The construction of martingales with given marginals has significance in finan- 
cial modelling (see for example [Q, |2|, [3j and |S]). In 0], the authors investigate 
the existence of an alternative model for which the Black-Scholes formula holds 
true. The case of the Bachelier formula raises the question of the existence of a 
non-Gaussian martingale with Gaussian marginals. 

Our solution to this problem is based on the elementary observation that for 
Y and £ independent standard Gaussian random variables, the distribution of 
Z = ^/rY + y/l — rt; is a standard Gaussian random variable for any value of 
r € [0, 1]. This allows us to randomize r and construct a family of Markovian 
martingales with Gaussian marginals. 

The question of constructing (Markovian) martingales with given marginals 
has seen considerable interest in recent years, mostly initiated by the paper of 
Madan and Yor (0). 

In [Jj], the authors give three different approaches: a continuous martingale, 
a time-changed Brownian motion, and a cosntruction that uses Azema-Yor's 
solution to the Skorokhod embedding problem. These constructions are applied 
to a number of special cases including the case with Gaussian marginals. But 
out of the three constructions, only the Skorokhod embedding approach yields a 
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non-Gaussian martingale. The other two reduce to a construction of a Brownian 
motion. The continuous martingale approach looks for a process of the type 



X t = [ a{X Sl s)dW s , 
Jo 



for which the marginal densities g(x, t) are iV(0, t). Writing the forward equation 
for these densities, it follows that a 2 = 1 (see |B]), and X t is a Brownian motion. 
In the time change approach 

X t = B Lt 

where L is an increasing process - in fact L is assumed to be an increasing 
Markov process with inhomogeneous independent increments, independent of 
the Brownian motion B. The assumption of Gaussian marginals implies 

and it follows that L t — t. But the Skorokhod embedding approach of Madan 
and Yor yields a discontinuous and time-inhomogeneous Markov process. 

Our approach is different to the above. It uses basic principles, and has the 
advantage of producing an entire family of processes indexed by an (infinite) 
family of subordinators. The construction produces a family discontinuous and 
time-inhomogeneous Markov processes. We obtain the quadratic variation of 
these processes and infinitesimal generators in some cases. The richness of the 
family has the potential to allow for the imposition of specifications other than 
the marginal distributions. 

Note that our method can be extended to include other types of marginal 
distributions, but for clarity of presentation we choose to focus solely on the 
Gaussian case. 

Finally, all existing approaches yield discontinuous processes (barring the 
Brownian motion itself), and the question of the existence of a non-Gaussian 
continuous martingale with Gaussian marginals remains open. 

2. A Family of non-Gaussian Martingales with Gaussian Marginals 

In this section we construct a (non-Gaussian) Markov martingale X t the marginals 
of which are Gaussian with mean zero and variance t. The existence of such 
process is guaranteed by a Theorem of Kellerer (see [H] and |BJ) that only re- 
quires the targeted marginal densities, g(x, t), be increasing in the convex order 
(E[/(X{)] > E[/(X S )] for s < t and / convex), and have means that do not 
depend on t. 

As eluded to in the introduction, the main idea of the proposed construction 
is the fact that, for any triple (R,Y,£) of independent random variables such 
that R takes values in [0,1], £ is standard Gaussian and Y is Gaussian with 
mean zero and variance a 2 , the random variable Z = a(vRY + a\/\ — R£) is 
Gaussian with mean zero and variance a 2 a 2 . However, the unconditional joint 
distribution of (Y, Z) is not bivariate Gaussian, as can be verified by calculating 
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the fourth conditional moment of Z given Y — 0. In fact, (Y, Z) is a bivariate 
Gaussian pair if and only if R is non-random. The martingale property of the 
two-step process (Y, Z) holds if and only if 



Y = E[Z\Y] = E[a(VRY + aVT^R£)\Y] = aE[\/R]Y, 



in other words, 



e[Vr] = 1 - 



(1) 



Furthermore, the conditional distribution of Z given Y is 



F: 



'z\Y=y{dz) = P[R= l]e ay {dz)+E <£ (oV%,aV(l - R),z\ 1 



L_R<1 



dz. 



where e x is the Dirac measure at x and 4>{p, a 1 , •) denotes the density of the 
Gaussian distribution with mean /i and variance a 2 . 

This construction of a two-step process can be extended to that of a con- 
tinuous time Markov process. Indeed, let i? Sj t be a family of random variables 
indexed by < s < t. We assume that R st takes values in [0, 1], has distribution 
that depends on (s,i) only through ■wt/s, and has moment of order 1/2 equal 
to \J s/t. We denote the distribution of R Sy t by G rrr{dr). We shall also need 



t/s K 

a family £ S] j of standard Gaussian random variables. 

The process Xt is constructed as a Markov process with the following almost 
sure representation of X t in terms of X Sl s < t, 



X t 




(2) 



We assume given the usual set-up of a probability space endowed with a filtration 
Tt to which X t is adapted. In the representation J2J, R s ,t and £ Sjt are assumed 
to be independent of each other, JF t -measurable and independent of T s . 

In the sequel, we shall often write a for y/s, a for y/t/s, t for yju/t and 
whenever independence between the variables involved need not be emphasized, 
R a for R St t and R T for Rt, u - 

Definition 1 The family {G a ) a >i is a log- convolution semi-group if the the 
distribution of the product of any two independent random variables with distri- 
butions G a and G T , is G aT . 

Define, for a > 1 and R a distributed as G a , U a = —\nR a , and, for p > 0, 
V p = U e p. If K p denotes the distribution of V p , then (G a )a>i is a log-convolution 
semi-group if and only if (K p ) p >q is a convolution semi-group: 



Proposition 1 Define, P St t(oc,dy) as, 

PoA x ' d v) z 



K p * K q = K p+q . 



2irVt 



exp 



2t 



dy, 



(3) 
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and for s > 0, 
P S;t (x,dy) 

= j(o-)eax(dy) + 



(4) 



1 



<*1> 1 -%tP^ )G a {dr)dy 



[0 ,i) V^\/iVT~- ~ "^ V 2i (! - r 
= j(cr)sax(dy) +E 0(CTv^"a;,a 2 CT 2 (l - R a ),y) 1r ct <i dj/, 

where R a is distributed as G a and 7(c) = G CT ({1}). 

J/ (Go- )<7>i is a log-convolution semi-group then for any u > t > s > and 
any £, 

P S) t(a:, dy)P ttU (y, dz) = P s , u (x, dz) (5) 



and, for any u > t > 0, 



iV(0,dy)Pt lU (y,dz) = P , u (0,dz). 



(6) 



Proof: We prove JSJ and JBJ) by showing that the almost sure formulation J3J 
is consistent. 



X u = r I y/Rt~uX t + OOL^Jl - Rt, u £t,u 



= T ( V^^ I VR^tXs + ay/1 - R s ^s 



t\j\ - Rt, u £t,u) 



= O-T \/Rs,tRt,uX s 



(1 — Rs,t)Rt,uS,s,t + v 1 — Rt,u£t,u 



Now, letting i? Si „ = R s ,tRt,u and 



?s,« — / J-fls „<1 + 1 -R S ,„=1 S«,t H / tfl s „<l?t,«) 



-X"u = TO- ( yjRs, u X s + ay/1 - R s ,u^s,u 



we see that 



with R su distributed as G T {dr). Also, the unconditional distribution of £ S)U as 
well as its conditional distribution given R Sj t and Rt lU are standard Gaussian. 
This in turn implies that £ Sjtl is independent of R s ,u- O 

Proposition 2 (Levy-Khinchin Theorem) Assume that the family (G a ) a >i 
is a log- convolution semi-group and let {R a )a>\ be independent random vari- 



ables with laws (G tT ) cr >i. Let L a (X) 



X\nR a 



E [(i? cr ) A ] be the moment 



generating function of the (positive) random variable U a = — In R a . 
For any a > 1, U a = — lni? CT is infinitely divisible, and 



In MA) 



/3A 



1-e 



-\x\ 



v{dx) 



lner 



(7) 
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/■oo 

where the Levy measure v(dx) satisfies ^({0}) = and I (1 A x)v(dx) < oo. 

Jo 

In what follows, we denote by ip the Laplace exponent of the log-convolution 

semi-group {G„) a >\'. 

/•OO 

ip(X) =/3\+ (1 - e- Xx ) v{dx). 

Jo 

The following theorem follows form the above and the Chapman-Kolmogorov 
existence result. 

Theorem 3 Assume that the family (G CT ) cr >i is a log- convolution semi-group 
with Laplace exponent 

/"OO 

ip(X) =/3X+ {l-e- Xx )v{dx). 

Jo 

If ip (1/2) = 1, then there exists a Markov martingale X t starting at zero with 
transition probabilities P Sy t(x 1 dy) given by J3J) and @) the marginal distributions 
of which are Gaussian with mean zero and variance t. 

3. Path properties 

Theorem 4 The process X t is continuous in probability: 
Vc > 0, limPNXt - X s \ > c] = 0. 

Proof: Using Lemma El below, we write, 

P[\X t - X,\ > c] < \v[(.Xt - X s f] = \[t- t 1 - 5 * 5 + t'-'s 5 - s] = t —^. 



□ 



Lemma 5 Let 5 = ip(l)/2 so that L a (l) = <j- 2S . Then 

E[(X t - X S ) 2 \X S ] = t - t x "V + t^s-^Xl - Xl 

Proof: Using representation J5J, we see that 

E[(X t -X s ) 2 \X s ] = V[V[(X t ~X s ) 2 \X Sl R a ]\X s ] 



,2 



aVE [1 - R a ] + E (ay/R a - 1 

oV (1 - L CT (1)) + (cr 2 E[i? CT ] - 1) X 2 
a 2 a 2 (1 - L CT (1)) + (a 2 L CT (l) - l) X 2 
a V - a 2 a 2 - 2 " + a 2 - 2S X 2 - X 2 



X 2 



D 
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Theorem 6 The (predictable) quadratic variation of X t is 



(X,X) t = St +(1-5) 



* XI 



ds, 



where 5 — ip(l)/2. Furthermore, it can be obtained as a limit in probability, 



n-l 



(X,X) t =Plim^E[(X tfc+1 - X tk f\X tk ] 

k=0 

where to < t\ < . . . < t n is a subdivision of [0, i\. 

Proof: First note that X t is a square integrable martingale on any finite interval 
[0,T]. In fact supE[X t 2 ] = T. Also, 

t<T 

E[X?\F S ] = E [t(l - R a ) + a 2 R a X 2 s \X s ] = t{\ - L a (l)) + a 2 L a (l)X 2 . 
Since L a (l) = a-- 25 = s 5 t- s , we find 

V[X 2 \T S ] = t - t^s 5 + t^s-^X 2 . 
It follows that 



E 



r X 2 
h-S) / -±du 
Jo u 

= (1 - 5) [ ^-du + (I - 5) [ (1 - u~V + u- s s- 1+5 X 2 )du 

JO u Js 

= (1-5) [ ^-du + (I - 5)(t - s) - s s (l - s _1 .X*)(t 
Jo u 



l v 2\(A-& _ s l-<5~) 



(1-5) 



s Y 2 



du + (1 - 5)(t - s) - s d t l - d + s + e- d s- 1+d X^ - X 



1-8 „ -1+5 y 2 



and 



E 



X 2 -5t-(l-5) 



t Y 2 

u 



l -du 



Ts 



s Y 2 



= t - t l - d s d + t l - d s- i+d X 2 -St -(1-5) / -^du -(l-5)(t- s) 

Jo u 

fl V-« _ s - t l - s s- l+5 X 2 s + X 2 

f s X 2 
= t-5t-(l-5) —±du-(l-5)(t-s)-s + X 2 
Jo u 

= x 2 - 5s -(1-5) r 

Jo 



X 



du. 



D 
The next result states that the only continuous process that can be con- 
structed in the the way described in Section |2] is the Brownian motion. 
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Theorem 7 If R a is not degenerate (R a ^ a~ ), X t is not continuous. 

Proof: We proceed by contradiction and assume that X t is continuous. Ito's 
formula for e lXXt gives 



s iXX * = 1 + M f - - ^-' xv 



e iXX °d(X,X) s , 



(8) 



where M t = / iXe l B dX s is a true martingale. In fact 



E[\(M,M) t 



- I X 2 e axx > 



d(X,X) ! 



- I X 2 e i2XX ° ( 5ds + (1 - S)'-^ds 



s 



E[Xl 



-ds 



= E 

= E 

< 6X 2 t + (l-6)X 2 

= SX 2 t + (1 - S)X 2 t 

= X 2 t, 



since X s is N(0, s) and ELY 2 ] = s. 

Taking expectations in ©, we obtain that 8(\t) = B[e iXXt ] = e~ A2 ' /2 must 
satisfy 



0{\,t) 



A 2 

1 

2 

A 2 



5 f 6(X, s)ds + (1-6) [ E [X 2 e lXX °] ds 
Jo Jo 

/•* d 2 e 

o d# 



,) I 9(X,s)ds-(l-S) f ^(X,s)ds 



Differentiating in t, we get that 0(X,t) must satisfy 



X 2 X 2 



SS(\,t)- 1(1 -s) w (\,t) 



that is, 



-Y = -y[*-(i-W 2 *-i) 



This, of course, can only occur if 6 — 1, which corresponds to Z CT (1) = cr 2 and 
i? CT being non-random equal to o~ 2 . □ 



4. Explicit Constructions 

Before we engage in the explicit construction of the processes outlined in the pre- 
vious sections, let us observe that these fall into one of two subclasses according 
to whether or not G a ({l}) is nil, uniformly in a. 
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Indeed, 

7 (cr) = G CT ({1}) = limL (7 (A) = lim cxp (-tp(X)lna) 

A"|"oc Afco 



and 



*y(<j) = <& lim ip(X) = +00. 

A| oc 



4. J. The Case 7(0-) > 

In this section we apply our construction to the case where 7(c) = G (T ({1}) > 0. 
The processes thus obtained are piecewise deterministic pure jump process in the 
sense that between any two consecutive jumps, the process behaves according 
to a deterministic function. Examples of such processes include the case where 
G„ is an inverse log-Poisson distribution. 

The interpretation of these processes as piecewise deterministic pure jump 
processes requires the computation of the infinitesimal generator. 

Proposition 8 LetG a be a log- convolution semi-group for which 7(17) = G a ({l}) > 
0, 7 is differentiable at 1 and lim -0(A) = 0. Then the infinitesimal generator of 

AJ.0 

X t on the set of Cq -functions is given by 
Aof(x) = 2^"^) and f° r s > °> 



[f{x + z)- f(x)} I 0((VF - l)x, 8(1 - r), z)G(dr)dz, 

[0,1) 



A s 


m 


= —fix) 

2s J V ' 




+- 


^/./ 


where 







G(dr) = lim 



G a {dr n[0,l)) 



.11 G ff ([0,l)) 
is a probability measure on [0, 1), and the limit is understood in the weak sense. 

Thus the process X starts off as a Browwnian motion and, when in x at time 
s, drifts at the rate of x/(2s), and jumps at the rate of — 7'(l)/(2s). The size 
of the jump from x has density /,„ ... <j>((y/r — l)x, s(l — r), z)G(dr), the mean 
of which is J, Q j\(\/r — l)G(dr)x. In other words, while in positive territory, X t 
continuously drifts upwards and has jumps that tend to be negative. In negative 
region, the reverse occurs; X t drifts downwards and has (on average) positive 
jumps. 

Proof: First note that the conditional moment generating function of U a 
given U a > is 

r* m _ L a {\)- 1 {(j) 

1-7(0-) 



and converges to 
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lim^(A) = l + ^. 
ct|i 7'(1) 



By the (Laplace) continuity theorem, if liniAxo "0(A) = then there is exists a 
probability measure on [0, 1), G{dr), such that 



G{dr) = lim 



G ff (drn[0,l)) 



<rii G CT ([0,1)) 



Next, 



f-s 



(E[/(X t )|X a = s] - /(»)) 



f(ax)j(a) - f(x) 



a 2 -I 



1 



+ ^T^ / /(») / ^(trVFar,t(l-r),y)G <r (dr)dy 
o" - 1 J JfO.l) 



[04) 

f(ax)j(a) - f(x) 



a 2 -I 
l- 7 (cr) 



(T 2 -! 



[o,i) 1 - 7W 



Letting ct decrease to 1, we see that 
AJ(x) 



1 



lim ■ 

t{a t — S 



(E[f(X t )\X s =x]-f(x)) 



xf'(x)+y(l)f(x) 7 '(1) 



f(y) I <f>{Vrx,s(l-r),y)G(dr)dy 

[0,1) 



= £/'(a0 + -^ / [/(y) - /(*)] / 0(VFs, s(l - r), y)G(dr)dy 



2s" 

, -Y(i) 

2s 



2s 



[0,1) 



[f(x + z)- f{x)} / 0((Vr - l)a?, s(l - r), ^)G(dr)dz. 



[0,1) 



D 
Note that the domain of A s can be extended to include functions that do not 
vanish at infinity, such as f(x) = x 2 . Indeed by Theorem^ g s (x) = 8+ (1 — #) — 
solves the martingale problem for f(x) — x 2 . 

The next proposition immediately follows from the observation that the pro- 
cess X does not jump between times s and t if and only if X u — «/^X s for 

u e (s,t). 
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Proposition 9 Let T s denote the first jump time after s > 0. Then, for any 
t> s, 

¥[T S >t]= 7 (a), 

where as before, a — yt/s. 



4-2. The Poisson Case: 7(0-) = cr 



1 



In this case f3 = 0, v(dx) = cSi(dx) with c = r-rr-, and ^p(X) = c(l — e ). 

1 — e i / 2 
In other words U a = — In i? CT has a Poisson distribution with mean c In cr. 

The assumptions of Proposition [S] are clearly satisfied with 7(c) = o~ c '. 

7'(1) = — c, limL*(A) = e~ A and G(dr) = e e -i(dr), so that X t has infinitesimal 

o-J.1 

generator 



^./(s) = ^/'(*) + ^ / [/(« + «) - /(«)]^(-«/c, s(l - e- 1 ), z)dz. 



It jumps at the rate of — with a size distributed as a Gaussian random variable 

x ~ 

with mean and variance s(l — e _1 ). The graph below shows a simulation of 

c 
a path of such a process. 




Furthermore, the law of the first jump time after s is given by 

,c/2 



»[T fl > t] = 7 (a) 



t c/2' 



In other words, T s is Pareto distributed (with location parameter s and scale 
parameter c/2 ~ 1.27). In particular, 



E[rj = 



c-2 



and E[T S 2 ] = 00. 
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4.3. The Case 7(0-) = 

When 7(c) = 0, we are only able to compute the infinitesimal generator for 
functions of a specific type. Examples of such functions include polynomials. 

Proposition 10 Assume that (3 — so that 



V(A) 



(1 



-\x\ 



v{dx) 



Let f be a C -function with the following property: there exist a function Nf 
and a (signed) finite measure Mf such that 

f (ae^x + ViVT^z) = *,(,) r e--M /(s , *, z , dX), u > 0, 



and 



Then, for any s > 0, 

A,m = f f{x) 

Is 



limiV/(a) = 1. 

o-J.1 



2.s 
Proof: Let 



1 /" /'+00 

+T" / [/(* + y) - fix)] / H(^ /2 - l)x, s(l - e- u ),y)v{du)dy. 



C„f{u) = C„f{s, x, z, u) = f (ae- u/2 x + ViVl - e~ u z\ . 

Then, since 7(c) = 0, U a is almost surely strictly positive and, 
1 



t-s 



(E[f(X t )\X s =x)-f(x)) 



1 1 

s a 2 — 
1 1 

1 1 



- T — - / (E f (<re- u ° /2 x + ^\-e- v °z\ - f{x)\ <P{z)dz 
(E [C*f{V a )\ - C a f(0))<f>(z)dz + (f(ax) - f(x)) 



S a 2 - 1 
f 1 



E 



N f (a) {e-^ - 1) M f (dX) 



(z)dz + (f(ax) - f(x)) 



- I N f (a) / / (e-^ 1 ™ l) M f (dX)<P(z)dz + (f(ax) f(x)) 



so ■■ - ± ( J ./„ 

1 1 f N f (o-)lno 



s o + 1 I o — 1 



'^'-'-MWMdz+f™-™ 



In o 



o-l 



Taking the limit as o J, 1 (that is t J. s), we get 

x 1 C f°° 

A s .f(x) = -f\x) --J J iP(X)M f (s, x, z, dX)cj>{z)dz. 
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Jo 



/>oc 

/ (l - e- Aw ) i/(dw) 
Jo 



Mf(s, x, z, d\)(p(z)dz 



(l-e- XuJ )M f (s,x,z,d\) 



v{duj)(j)(z)dz 



f [e- u/2 x + y/sVl - e~"z\ - f(x) v{du)<j){z)dz 



and the proof is completed by a change of variables in z. 
Lemma 11 Let f{x) = x n , then 

f (ae- u/2 x + \/iVl - e~ u z) = a n f e _Au M / (s, x, z, dX) 

where 

n 

M f (s,x,z,dX) = J2Yl 



D 



,' 



k\j\(n — k — j)\ 



(-iyx k s ^- k ^ 2 z n - k {e k/2 * m] ){dX) 



fe=0 j=0 

and mj(dX) is the j -order convolution of the probability measure 

m(dX) = -^^ T{n -^ /2) s n (dX). 

2\/ir z — ' n\ 



v 71=1 

Proof: First, write the Taylor series of the (analytic on (0, 1)) function 1 

l£Sr(n-l/2) „ 



VI - x, 

2 ^ n!r(l/ 
It immediately follows that, 



ti " !r (V2) 



i-vT 



l^r(n-l/2) 



£ 



1 n!r(l/2) 



, ' --== / e" A "m(dA), 
o 



+oo 



1 ~r<-*j -p/ -i /q\ 

where m{dX) = r _ > j e„(dA) is a probability measure. Now, 



2^^ n! 

v n— 1 

/ fcre-" /2 a; + \/t\/l - e- u z) 
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/V ( n \e- ku ' 2 X k S ^- k ^ 2 (l - e -u)(n-k)/2 z n-k 

fe=0 ^ k ' 
n 

n 

- B EE 



fe=0 

fe=0 
n n — k 



e -ku/2 x k s (n-k)/2 



i — k 



n — k 



k=0 j=0 



k\j\(n — k — j) 



1 - (l-Vl-e- 



The proof is ended by observing that 



-ku/2 ^ _ y/iZ^y = r e ~ Xu (e k/2 * mj ) (dX) 



D 



The following theorem is now proven. 



Theorem 12 Assume that 3 — 0. For any polynomial f and any s > 0, 



A s f(x) = -f'(x) 
2s 

+ Ys J [I{X + V) ~ f{x)] J ^ e ^ /2 - ^ s{1 - e-n,y>(du)dy.(9) 



4-4- The Gamma Case: 7(c) = 



Here 3 = 0, u(dx) = ax 1 e bx dx with a = — -, =-r- and %b(\) = a In ( 1 + — 

ln(l + ^) ^ V V & 

that is C/ CT has a gamma distribution with density 



h a (u) 



b" 



lnr 



r(alno") 



f alna-l e -6« U>Q 



and _R CT has an inverse log-gamma distribution with density 



La In a 



r(r) 



r(alncr) 



(_ lnr )aln ff -l r (,-l ; <r<l. 
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In this case it is possible to compute the generator for a much wider class of 
functions. 

Proposition 13 Let G a be the log-convolution semi-group of the inverse log- 
gamma distributions. Then \Q holds for any bounded function with bounded first 
derivative. 

Proof: In the proof of Proposition ^J we write that 

(E[f(X t )\X s = x] - /Or)) = \^— [ f (E \C a f{U a )] - f(x))<f>(z)dz. 



t-s 



Denote by 9(u) the quantity e u / 2 x+i/i\/l — e~ u z. Then, inserting E [Cif(U a 
E[/(0([/ CT ))]weget, 



t-s 



Since 



(E[f(X t )\X s =x]-f(x)) 



1 1 



SO-+1 



,[C a f{U„)]-nCif{U a )] , E[d/([/ CT )]-/(.x) 



o--l 



a - 1 



(10) 

4>{z)dz. 



c a f{u a ) - Ci/om f(o-e(u a )) - f(0(u a )) 



= eiu^fM, 



£7-1 (7-1 

for some r\ a between 9(U a ) and aO{U a ). 9 and /' being bounded, we obtain that 

\[c a f{u a )]-ncif(u a )] 



lim 



C7-1 



(z)dz = xf'(x). 



To compute the limit of the second term in (|10|) . we use Lemma IT3I below . which 
shows that 

-11 / C7-1 ^ W 
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, , f {e- u ' 2 x + Jsy/l - e- u z) - f(x) , 



o 



[/ (x + y) - /(a;)] / <j>{x{e~ u / 2 - 1), s(l - e~ u ), y)a dudy. 

Jo u 

D 
Note that since i/((0,oo)) = +oo, f °° <j)(x(e~ u / 2 — l),s(l — e~ u ),y)i>(du)du 
cannot be re-scaled to produce a density for the jumps of the process. 

Lemma 14 Let V p have a gamma distribution with density: 

KP 

Let g be such that g(0) = and g(v)/v is bounded. Then 

lim-E[g(V p )} = r ?^±e- bv dv. 
PiO p Jo v 



Proof: First observe that 



-E[g(V p )] = U 
p b 



g(v P+ i) 



V P+ i 
taking the limit as p J, 0, we obtain by dominated convergence 



limiEb(y p )] = il 



9(Vi) 



Vi 



g(v) - bv . 

e dv. 



a 
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